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k*" , Abstract 

We discuss the possibility of extracting \ V c b\ from the future measurement at LHC of the B c — ► J ftp fi + v 
decay rate followed by the leptonic decay J/?/' — > /J + [i~ ■ We examine in detail the influence of the mea- 
surement of the muons' momenta and vertex reconstruction in the final uncertainty of \V c t\, as well as 
the required theoretical inputs at both the production and weak decay level, and their reliability in this 
regard. We stress the relevance of a rigorous description of the fragmentation contributions to inclusive 
production at large transverse momentum of heavy quarkonia in pp collisions. 
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a: 

^ • 1 Introduction 



The interplay between theory and experiment, inherent to all branches of science, be- 
comes crucial in high-energy physics at present. Indeed, despite any (always welcome) 
unexpected discovery, the required investment of money and time demands well-defined 
aims and physical objectives from the theoretical side long time before the start of real 
data-taking. This is particularly apparent at LHC experiments, in one of which (ATLAS) 
one of us (M.A.S.L.) is currently involved. On the other hand, it is a task of theorists to 
survey all those theoretical topics compatible with the ultimate goals of an experiment, 
keeping connection with a "realistic" experimental point of view. This is, in fact, one of 
the underlying motivations of the present work. 

In particular, it has been recently pointed out the feasibility of the observation of B c 
mesons in the ATLAS experiment With respect to their production rate, setting the 
bb cross section equal to 500 fib and assuming the fragmentation probability of a b quark 
into a B c meson of the order of 10~ 3 0, the expected number of B c mesons at "low" 
luminosity [] is about 10 10 per year (10 7 s) corresponding to an integrated luminosity of 
10 fb- 1 . ' 



* Research partially supported by CICYT under grant AEN 93-0234 
' galdon@evalvx.ific.uv.es 

* mas@ evalvx .ific.uv.es 

x It is foreseen that during the first years, LHC will run at C ~ lxlO 33 cm -2 s _1 . Later on, the luminosity will increase 
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In this paper we shall argue that in LHC experiments the decay channel 

B c -> J/^j ^+ V p (1) 

followed by the leptonic decay of the J/ip into a pair of oppositely charged muons, could 
provide a reliable extraction of the mixing matrix element \V c b\. The expected semileptonic 
branching ratio ~ 2% for channel (1)0 together with the leptonic branching ratio for the 
J /if) into two muons ~ 6% @, yields an overall branching fraction of order 10 -3 . 

We next combine the last BF with the acceptance Q for the muon trigger ~ 15% (where 
the three possible triggering muons per decay have been taken into account.) Besides, the 
detection of the two remaining particles (i.e. the non-triggering muons satisfying some p± 
and rj cuts ||) amounts to an acceptance for signal events of ~ 12%. Finally, we assume 
an identification efficiency for each muon of 0.8 H yielding a combined value of 0.8 3 ~ 0.5. 

Thereby, the expected number of useful B c decays (] could reach several 10 5 per year at 
low luminosity. 

Such large statistics together with the foreseen precision in the muon momentum mea- 
surement should permit the experimental access to the kinematic region near zero-recoil 
of charmonium (with respect to the B c meson) with a good energy/momentum resolution. 
Moreover, very stringent cuts can be put on events and hence, expectedly, background 
could be almost entirely removed from the event sample permitting a very precise deter- 
mination of the B c lifetime. We explore all these issues in the following section. 



2 Extraction of \V c b\ 

In recent times, an attempt to derive the value of the Cabibbo-Kobayashi-Maskawa matrix 
element \V c b\ from B^D*iV in a (more or less) model- independent way has been proposed 
within the framework of the heavy quark effective theory (HQET)||. This method 
basically relies on the existence of a universal form factor, the so-called Isgur-Wise function 
£(w), depending onto = frt>2, variable representing the Lorentz factor of the final hadron 
in the initial hadron's frame. (Hereafter subindices 1, 2 will denote initial, final hadronic 
quantities.) 

Near zero recoil 



1 dT(B^D* en) Gl . x9 o 9 ? 9/ . , l2 

— = 4# (mi ~ m2) ml Va m lVcbl (2) 

where t)a denotes a perturbative coefficient |5|. On the other hand, the i(w) form factor 
absorbs those corrections arising from inverse powers of the heavy quark masses, thus 

2 Hereafter, we mainly rely on technical details contained in the Technical Proposal of the ATLAS Collaboration ||. Those 
aspects directly related to the generation of the B c signal and its background, the trigger and reconstruction efficiencies are left 
to a separate publication Let us only mention that first-level trigger muons are required to have a transverse momentum 
higher than 6 GeV/c and pseudorapidity |?j| less than 2.2. The other two muons are required to pass the cuts p± > 3 GeV/c 
and \rj\ < 2.5 

3 Semi-electronic B c ^ J/ipe + v decays, followed by the J/ip decay into fi + [i~ can be considered as well. In addition, if the 
muon from the semi-muonic decay (1) gives the trigger, the event can be selected with the mode J/tp^ e + e~. The trigger 
acceptance for the former is ~ 12% whereas for the latter is ~ 3% [H| 
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removing its universal character. In the infinite quark mass limit, £ coincides with the 
Isgur-Wise function. 

Experimentally, one can obtain the product £(u>)-|V^,| from the measured differential 
rate (l/Te)dBr /dq 2 , where q = P1—P2 is the four-momentum transfer to the leptonic sys- 
tem (see figure 1). As developed by Neubert [0], the zero recoil point is especially suitable 
for the extraction of |T4&| in B^D*£v decays. This is because heavy-quark flavour sym- 
metry and Luke's theorem determine the normalization of £(1) up to second-order power 
corrections: £(1) = 1 + 5\/ m 2 where the last coefficient stands for the long-distance cor- 
rection ||. The theoretical control of these uncertainties permits a reliable determination 
of \Vcb\ at this kinematic point where, however, there are no data due to vanishing phase 
space. An extrapolation to this endpoint is thus required and the final precision depends 
both on the experimental accuracy, background removal and theoretical corrections as 
well. 

In order to apply a similar approach to the semileptonic channel (1) of B c mesons we 
will characterize the kinematic state of hadrons by means of their four-velocity instead 
of four-momenta as well. This will make sense in our later theoretical approach as we 
are dealing with weakly bound states where the constituent heavy quarks are not too far 
offshell, as pointed out elsewhere [IH] [TT| . 



For clarity of presentation let us make below some remarks to be developed one by one 
in the following: 

a) The accessible w range in B c — > J/ip^u decays (see figure 2) is limited to the interval 
(1,1.26), whereas in B-^D*lV decays the range is significantly wider (1,1.5). The 
full kinematic reconstruction of the decay is required in order to measure w for each 
sampled event. 

b) The precision in the measurement of the muon's momentum in 5-physics at LHC 
experiments is of order 2% || or less |§ , representing an uncertainty of about 2% or 
less in the reconstructed momentum of the J/ip. 

c) Background from processes like B^J/ipK^ must be taken into account. 

d) Contamination from B c ^Xfi + i> decays, where X = ip(2S), Xc(l-P) stands for a higher 
charmonium resonance, will be addressed as well. 

e) We shall consider combinatorial background coming from bb pairs subsequently de- 
caying into muons. 

f) A rigorous fragmentation framework recently developed by Bodwin, Braaten and 
Lepage |nj could provide a reliable evaluation of the (high-pjj production cross 
section of heavy quarkonia in pp collisions at LHC energy. 

g) Precise predictions for the decay rate in the vicinity of zero-recoil could be obtained 



in the framework of lattice QCD |T5[ or even potential models |16| of hadrons. 
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2.1 Signal 



In this section we develop items a) and b): 

a) In order to get the w differential distribution of the B c semileptonic decays, their 
complete kinematic reconstruction (event by event) is required. This can be achieved 
despite the unknown energy of the initial (B c ) particle (whose mass is assumed to have 
been previously measured from the B c ^J/^n decay at the LHC IJ, or even earlier at the 
Fermilab Tevatron) and the unknown momentum of the final v, if the initial direction of 
the former can be experimentally determined. Indeed, this should be the case by means 
of the vertexing capability of high performance inner detectors at LHC experiments. 

In fact, a precise secondary vertex resolution for B physics has been established as an 
experimental goal in the ATLAS experiment || 0. This fact together with the accurate 
reconstruction of the primary vertex at low luminosity should permit the measurement of 
the Bc's flight-path . 

Let us write w as 

w = 7i7 2 (1 — \v 2 \cos6) (3) 

where 7« = (1 — |-Uj| 2 ) -1 / 2 , stands for the modulus of the three- velocity of particle-i and 
6 denotes the angle between the J/ip and B c (reconstructed) trajectories. The influence of 
the geometric uncertainty on w through the cosine in the above quantity near zero-recoil 
(7! ~ 72 = 7) can be inferred from the following expression 

5w ~ \v\ 7 y/2 (w - 1) 56 = ^2 (w - 1) 56 (4) 

where 56 = a/(n+ 1)j\v\t and 56 = cr/(n + l)cr. We shall take o pa 100 /mi as 
a reasonable and safe estimate for the combined primary+secondary vertex resolution. 
(The vector joining the production and decay vertices in space could be measured with a 
precision of (of the order of) 100 fim.) The B c average decay length has been expressed 
in (4) as n + 1 "lifetime-units" where n denotes the cut on the proper decay length, i.e. 
we consider only decays with proper times larger than nr. We shall set t pa 0.75 ps as 
the expected -B c 's lifetime. This is halfway between potential model predictions (0.5 ps) 
and QCD sum rules expectations (1 ps) |T7[ . 

It is quite remarkable that the accuracy of w becomes independent of kinematic factors 
since a more precise measurement of the angle due to a longer decay length because a 
larger Lorentz factor, is compensated by a greater influence on 5w. On the other hand, 
notice that the uncertainty 5w increases with w. This behaviour may be qualitatively 
understood since a lower recoil implies a smaller 6 angle (6 = in the limit w = 1), thus 
a smaller uncertainty through cos 6. 

Tentatively we shall set n = 7 (thus 56q ~ 56 mrad) amounting to a reduction factor 
10~ 3 of the original sample of triggered and reconstructed events. In fact, this choice 
represents a compromise: the requirement of longer B c flight-paths by imposing higher 

4 It is foreseen the installation of an additional precision layer at low radius during the low luminosity running 



4 



n values should certainly diminish the experimental uncertainty on 9 but would imply a 
further loss of statistics, thereby increasing statistical errors. 

b) With regard to the influence on w from the muon momentum resolution we obtain 

5w ~ 2 (w - 1) — (5) 

7 

Setting £7/7 ~ 2%, we notice that numerically its actual influence on w is smaller than 
that due to the uncertainty through the angle 9. 

2.2 Background 

We address below those aspects related to contamination of the signal (1) by either B 
decays or combinatorial background: 

c) Background from decays like B^J /if)K ± may fake the B c signal if the K is misiden- 
tified as a muon. This could happen due to punchthrough, for example, if the kaon decays 
to a muon in flight before it reaches the calorimeter. Let us remark, however, that such 
muon could be rejected from the stand-alone measurement of its momentum in the air-core 
toroid muon detector of ATLAS, once compared to the (larger) momentum of the parent 
kaon independently measured in the inner detector. 

In addition, there are some possible strategies to remove this kind of contamination by 
firstly requiring coplanarity among the direction of the parent meson, the reconstructed 
trajectory of the J /if) and the n track. In fact, this requirement amounts to no missing 
transverse momentum with respect to the direction of the decaying particle. Besides, the 
invariant mass of the final particles would be compatible with the B decay hypothesis. 

There is, however, another potential source of contamination from B—>J/if)K* decays 
occurring again at a comparable rate as the signal. The K* subsequent decay into 
and an undetected tt° would make useless the last invariant mass condition. Nevertheless, 
since the momentum of the in the K* rest frame is 291 MeV 0, a lower cut of about 
300 MeV on the component of the momentum perpendicular to the plane defined 
by the decaying particle and the J /if) should remove all that background with a good 
acceptance (~ 80%) for the signal ||. In fact, this cut is compatible with the precision on 
the transverse momentum due to the flight-path reconstruction. Likewise, other sources 
of background from channels like B— ^XJi^ 1 where X c denotes charmonium resonances 
subsequently decaying into J /if) and neutrals, should be rejected as well by means of this 
cut. 

d) Contamination coming from B c ^Xfi + u decays where X = if>(2S), x c (lP), followed 
by the resonance radiative decay into a J /if), would contaminate the signal (1). Indeed, 
one cannot get rid of this contamination by imposing to the semileptonic decays a missing 
mass constraint, contrary to B—>D*iv in a e + e~ collider. However, it may be possible 
to detect the photon in the electromagnetic calorimeter, with typical energy of few GeV, 
predominantly emitted within a certain cone along the direction of the J /if). In this 
way those events yielding an invariant mass compatible with any of the above-mentioned 
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resonances would be removed from the event sample. Nevertheless, before any serious 
simulation analysis about this possibility had been performed, we prefer to stay at more 
conservative grounds and consider provisionally the emitted 7 as useless. Therefore we 
must evaluate the degree of contamination of the signal sample. 

We tentatively adopt experimental data on B decays by making the reasonable as- 
sumption that the B c should decay semileptonically into higher resonances and J /if} with 
a relative ratio of order 2.7/7 M. Further, using crudely an "average" BF ~ 0.16 for 
the decay of the charmonium resonances into J/ ip and neutrals , we conclude that the 
expected contamination of the signal from all the above channels represents about 6% 
hence below the statistical fluctuations of the sampled signal events as we shall see. 

e) Combinatorial background is generated by the simultaneous coexistence of three 
muons in events. For instance this could happen if a J /if} is produced from a B whereas 
the other muon comes from a different b quark (thus from two spatially well separated 
vertices). This kind of contamination should be drastically suppressed by means of the 
seven-lifetimes cut imposed to the reconstructed secondary vertices. Furthermore, prompt 
J/^'s produced at the primary pp interaction should be practically entirely removed as a 
background source by requiring a combined cut on the \ 2 °f the fitted secondary vertex 
and the decay length |J. 

2.3 Theory inputs and determination of \V c b\ 

Next we examine in detail the possibility of extracting \V c b\ from experimental data. 

f) Let us start by recalling that 

dT(B c ^J/jj fj&u) = 1 dBr(B c ^J/jj fj,+u) 
dw t dw 



The differential branching ratio dBr/dw can be written as 

dBr(B c ^J/ip n+v) 1 dN(B c ^J/ip ^v) 

dw ~ N{B C ) ~chL ^ ' 

where dN/dw represents the number of semileptonic decays per unit of w for a certain 
integrated luminosity, once corrected by detection efficiency; N(B C ) is the corresponding 
total yield of B c 's produced for such integrated luminosity. (Notice that r could be 
determined accurately from the same collected sample of semileptonic B c decays.) 

Let us observe that the exclusive decay mode B c ^J/ipn to be detected at LHC J! 
cannot be employed for normalization purposes as far as its BF would not be (obviously 
yet) experimentally known. On the other hand, B c decay modes involving the c-quark 
(which in principle would permit the measurement of ratios of CKM matrix elements in 
combination with b decays) will be hardly observable among the huge background of a 
hadron collider. 

Therefore, in order to get N(B C ) it will be convenient to compare the B c production 
rate to the prompt if}' yield in pp collisions, for the same integrated luminosity. Therefore 
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a(pp^B c + X) AT(^-,/iV 



a(pp^' + X) ep-BFty'^n+fi-) 

where N(tp'—>-fi + fi~) stands for the number of prompt ip' experimentally observed through 
the muonic decay mode and e^> denotes its detection efficiency Contamination from weak 
decays of B mesons into ip' states should be efficiently removed by means of the vertex 
capability of the inner detector fj • 

Therefore we suggest normalizing the B c sample with the aid of the ip' yield through 
expression (8). The ip' state is preferable to the J/ip mainly because the former expectedly 
should not be fed down by higher charmonium states Instead, data released by Teva- 
tron shows that the J/ip indirect production through \ c intermediate states may easily 
overwhelm direct production. Let us also remark that the production of the ip' resonance 
in pp collisions is more than one order of magnitude larger than initially expected ]l{J 



At sufficiently large p± f], one reasonably expects that the main contribution to heavy 
quarkonia production comes from the splitting of a heavy quark or a gluon. Indeed, 
even though the fragmentation process is of higher order in a s than the "conventional" 



leading order diagrams pTJ, the former is enhanced by powers of p±/rriQ relative to the 



latter [[| | 22|| . Although the situation is still controversial in the literature, there are 



strong indications from very detailed calculations (see [23 ]24|] and references therein) that 
fragmentation indeed dominates for large enough p±. In the following, we first examine 
some color-singlet mechanisms, expected to contribute largely to the high-p^ inclusive 
production of prompt heavy quarkonia. 

According to the set of papers in 0, perturbative QCD can provide a reliable calcu- 
lation for the fragmentation function of a high-p^ parton into heavy quarkonia with only 
few input parameters: the charm and bottom masses, and the square of the radial wave 
function of heavy quarkonium at the origin |i?(0)| 2 . Note that the largest uncertainty for 
normalization purposes comes from the (third power of the) charm mass in the fragmen- 
tation function (see @ for explicit expressions). However, we are only interested in the 
ratio 

= a(pp^B c + X) 

a(pp^' + X) v ; 

where each cross section can be written as a convolution of the parton distribution func- 
tions of the colliding protons, the cross sections for the hard subprocesses leading to 



the fragmenting gluon or heavy quark and the respective fragmentation functions || | 25fl . 
Thus notice that the fragmentation functions D^ B , D g _+ Bc and D c ^>, D g ^> f26|j , which 



themselves are independent of the parton-level subprocesses, appear combined as a ratio. 
Therefore, those uncertainties coming from the common overall factor vn? c automatically 
cancel each other. Moreover, those uncertainties introduced by the parton distribution 
functions also should significantly diminish in r as well, especially when a cut on the 



5 See Ref. jjjj] for an alternative explanation of the observed ip' surplus found in Tevatron. However this mechanism would 
require an uncomfortable large radiative BF for higher resonances 

6 We find an average p± « 20 GeV/c for B c mesons simply passing the kinematics cuts on p± and |tj| of the decay muons |H| 
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transverse momentum of heavy quarkonia is required. Still one must evaluate 



I^B C (0)| 2 flQl 

In fact there is the possibility of expressing Eq. (10) according to the general factorization 
analysis of [13] in a more rigorous way as 

<0\O?°( 1 S )\0> 

fto = n U J- 

<0|Of( 3 ^)|0> v ; 

where are local four quark operators and the matrix elements can be evaluated from 
NRQCD. A similar expression holds for the P*. (Moreover, the denominator can be 
determined from the measured leptonic width of the ip' .) 

Recently, an additional color-octet fragmentation mechanism [^7J] has been suggested 
in order to reconcile the experimental results on inclusive ip' production with theoretical 
predictions. This mechanism assumes the creation from gluon fragmentation of a cc pair 
in a color-octet state, in analogy to \ c production [jnj. Then a new nonperturbative 
parameter is required, appearing in r combined as the dimensionless factor 

KKV) ~ |iV(0)| 2 " <0|Of(35 1 )|0> ( } 

where we remark that it can be expressed in terms of NRQCD matrix elements |14j \\ 

Furthermore, if we do not neglect the contribution to the total P c production of those 
orbital excitations like P-wave states, new fragmentation functions -D£->6c(P) come into play 



30fl . Accordingly, some new /tj factors will appear which can be conveniently expressed, 



for instance in analogy to Eqs. (11) and (12), as: 

\R Bc (0)\ 2 K ™ 2 <0|Of e ( 1 5 )|0> 



2™ 2 H lCbc) 1 <0|Of-(P J )|0> 



and 



_ Wg 8(fe) ^ i <o\oj^s 1 )\o> 
( c) l^ c (o)l 2 ~™ 2 <o|of*(%)|o> li4j 



where rh is the reduced mass of the b and c quarks. The long-distance Pi parameter [2S 



is related to the square of the derivative of the color-singlet wave function at the origin . 
Higher B c resonances like P-wave states might be further taken into account, introduc- 



ing new nonperturbative parameters [31] involving higher derivatives of the wave function 



at the origin, or the NRQCD matrix element analogues. 

In sum, the evaluation of the Ki parameters would allow one to complete the compu- 
tation of the ratio r and thereby to obtain from the experimental measurement of the if)' 



7 -^8(ip') can b e phenomenologically estimated from B or T decays as well [p3 J2!|. On the other hand, according to a 
nonrelativistic quark model Hg is related to a fictitious color-octet wave function at the origin 
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yield, the total number of B c events for an integrated luminosity with the aid of expression 
(8). 

We next address the hadronic transition between the initial B c and final J/ ip particles. 

g) From the theoretical side the hadronic transition B c —>-J/iIj involves two heavy-heavy 
systems. Although originally HQET was only applied to heavy-light bound states, there 
is a growing belief that it can be cautiously applied to the former [|K|] Indeed, 



spin symmetry should be still valid as a first order approach in the dynamics of heavy- 
heavy hadrons. This allows the introduction of the analogue of the Isgur-Wise function 
for transitions involving doubly heavy hadrons, denoted by 7/ 12 (t>ri>2). At the non-recoil 
point v i = v-i — v , we shall write [TU| 

< J/$\ A» \B C > = 2 7? 12 (1) y/mmZe? (15) 

where A 11 = 07^756 stands for the axial- vector current and e 2 represents the four- vector 
polarization of the J/i[). 

Following similar steps as in the B decay into D*, we write the analogue of equation 
(2) as 

lim^-vi— ■= - —-{m 1 -m 2 ) m 2 T] 12 {L) \ V cb \ [Lb) 

y/w 2 — 1 aw Atx 6 

It is of key importance from the theoretical side in our proposal, the existence of a single 
form factor in the above expression to be determined theoretically in a rigorous manner. 
Observe that this is valid only in a neighbourhood of w — 1, however. Deviations from 
the strict requirement of the spin symmetry lead to the appearance of new form factors 
in the hadronic matrix element |§. Even though they can be absorbed in a redefinition 
of 7712 (w) this does not mean that we could ignore them if \V c b\ has to be extracted away 
from zero recoil. 

Note that 7712(1) may be interpreted in first approximation as the overlap of the initial 
and final hadron wave functions [|10] Due to the intrinsic nonrelativistic nature of 

its heavy constituents, nonrelativistic QCD on the lattice or even refined potential based 
models of hadrons, could be applied in order to get 7712 and its slope at zero recoil in a 
reliable way. 

This fact in conjunction with the presumed accuracy of the J ftp momentum measure- 
ment suggests determining \V c b\ in a similar way as in the method proposed by Neubert 
through the B->D*£V decay. 

Observe that from an experimental point of view the D* kinematic reconstruction 
should suffer from larger uncertainties through the cascade decays: D*^D 71 followed 
by the D decay into final states with only charged particles. Conversely, the J /if) decay 
into fi + fi~ is much cleaner by all means. Of course, the environment is much dirtier in a 
hadron collider than in a e + e~ factory 0. Nevertheless, since the signal consists of three 
muons the B c decay can be disentangled from the huge hadronic background by means of 
adequate kinematics cuts and constraints 0, at least at low luminosity. 



3 One may conjecture about the possibility of producing B C B C pairs in a future fj, + jj, collider 
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In order to make a crude estimate of the slope of the analogous Isgur-Wise function r]i 2 
at zero recoil, we shall make use of a nonrelativistic approach. Specifically we employ the 
ISGW |M| model which should give sensible results in the case under examination. The 
overlap wave function then reads 



\ 2/3^ - 


3/2 


U 2 + & 2 - 


exp 



Tfl t m 



2 mi m 2 k 2 {(3l + (3l) 



(17) 



where t = q 2 , t m = {mi — m 2 ) 2 and (3i denotes a variational parameter corresponding 
to the i- meson obtained by adjusting the hadronic spectra; A; is a parameter introduced 
somewhat by hand, usually set equal to 0.7. 

Now let us introduce the slope parameter p\ 2 of the rji2(w) function near zero recoil as 

Vl2 ( w ) = Vl2 (l) { 1 - p\ 2 (w-1) + 0[(w-l) 2 } } (18) 

Since t m — t = 2mim 2 (w — 1), we find that p\ 2 can be expressed according to the ISGW 
model as 



in 



2 



* - w, (19) 

where (3\ 2 = k 2 (/3 2 + /3f)/2. For transitions between heavy-light mesons, (3\ = (3 2 = (3 
and putting k = 1 one quickly recovers the result of ||35|| : p 2 = m 2 /2/3 2 . (The relativistic 
correction due to the boosted four momentum of the light quark q should not affect p\ 2 
so much.) 

Setting the parameters in Eq. (19) equal to some typical numerical values for B c and 
J/ij) states, (3\ = 0.82 GeV, (3 2 = 0.66 GeV, m c = 1.5 GeV and letting k vary between 0.7 
and 1, we estimate a range 

p\ 2 « 2 - 4 (20) 

representing a higher value than p 2 (of order unity) for transitions between B and D* 
mesons [0], as otherwise expected |36| . (Of course, this range of p\ 2 values must be 
considered as merely indicative.) Note also that one expects the curvature of r]i 2 {w) at 
w = 1 for transitions between heavy- heavy hadrons to be larger than that of £(w ) between 
heavy-light ones. This implies the necessity of a closer access to the zero recoil point for 
the accuracy of the Taylor expansion of Eq. (18) as is indeed the case. 

For the purpose of illustration, we have depicted in figure 3 some expected typical 
points and error bars for the decay B c ^>J/ifj(—>p + p~) p + v for a definite prediction of the 
ratio of production cross sections r in Eq. (9). 

Under this condition, let us remark that the uncertainty on w (i.e the horizontal er- 
ror bars) is of systematic nature, mainly coming from the B c flight-path reconstruction. 
With regard to the vertical coordinate rji 2 (w)-\ Vd>\, we have assumed that its uncertainty 
is essentially statistical, besides the contamination by decays into higher charmonium 
resonances which are included as well in the vertical error bars. 
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3 Summary and discussion 



This work concerns B c mesons to be copiously produced at LHC next century. We have 
argued that it should be possible to reconstruct those semileptonic decays into a J /if) near 
zero recoil, followed by its consequent decay into a muon pair. The fact that the J/ip is a 
vector particle should lead to a mild falling off of the number of events in the kinematic 
region near zero recoil as in B^D*£u. (This is because of the s-wave contribution to these 
decays, absent for a final pseudoscalar hadron as in the B—>Dtv channel.) 

The signature of the decay B c ^>J/if)(-^fi + fi~) /i + z/, consisting of three muons coming 
from a common secondary vertex could be disentangled from the huge hadronic back- 
ground. The excellent accuracy in the measurement of the muons momenta would lead to 
a precise determination of the J/ ip momentum and energy. In addition, a good 3D vertex 
reconstruction capability is determinant for the full kinematic reconstruction of events. 
We have shown that, in fact, it is 58q = a/(n + l)cr one of the parameters that actually 
matters in the expected accuracy of the differential w distribution of sampled events. 

We also addressed some background sources, either from real B decays but faking B c 
ones, or from true B c semileptonic decays into higher charmonium resonances. Several 
strategies to remove all these contaminations were proposed. 

From the theoretical point of view, the heavy quark spin symmetry permits the in- 
troduction of the analogue of the Isgur-Wise function r/u for transitions between doubly 
heavy mesons. Deviations from the spin symmetry imply the appearance of new form 
factors in the hadronic matrix element which, however, do not contribute at zero recoil. 
Hence, at this kinematic point, only a single form factor (though including itself spin 
breaking effects) has to be determined for the B c -^J/ip transition. 

Experimental data would permit to obtain the product t]i2(w)-\V c b\ near zero recoil in 
a more or less accurate way according to the accuracy of the differential decay rate dT/dw 
itself. The measurement of the latter requires the normalization of the B c semileptonic 
events to the total number of the B c sample. (In other words the knowledge of the absolute 
branching fraction for the semileptonic decay is required.) To this end, one can rely on the 
fragmentation approach providing in a rigorous way the ratio of the yield of B c mesons 
relative to the prompt ip(2S)'s produced in pp collisions (see Eq. (8)). As far as we 
are interested in their relative production rate r, the uncertainty should be much less 
dependent on theoretical or phenomenological inputs (like heavy quark masses) than in 
each production rate separately. 

We have shown that a set of long-distance dimensionless factors (denoted by Ko,i,s) 
involving ratios of nonperturbative parameters like |i?(0)| 2 , Hi and H' 8 must be determined 
to make a reliable prediction of r. Even if their evaluation from NRQCD [|14j] is still in 
its infancy, lattice techniques should permit in a not too far future to go beyond model 
estimates, including relativistic corrections, in a systematic and reliable way. 

In addition, forthcoming e + e~ factories and Tevatron running at higher luminosity 
should improve considerably the accuracy of the measured BF of the B and T decay 
modes involving P-states of charmonium, improving the experimental information on such 
parameters. Let us also remark that the fragmentation hypothesis can be checked with 
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experimental data from pp or pp collisions, for instance through the transverse momentum 
dependence of heavy quarkonia, its polarization and association to jets of light hadrons. 

Once the differential rate dT/dw is achieved to be experimentally known, an extrap- 
olation to the w = 1 endpoint, either linear or quadratic (i.e. including the curvature) 
would yield r/i2(l)-|Vd,|. Consequently, the knowledge of 1/12(1) should provide \V c b\. Let us 
stress that this calculation and the evaluation of the fragmentation probability into heavy 
quarkonium, both involving heavy quarks, rely essentially on similar theoretical grounds 



and could be derived from first principles [14 



Of course, one may consider as well the possibility of choosing a kinematic point 
different from w = 1 to obtain \V c b\- This requires, however, the introduction of new 
independent form factors (i.e. entering in the hadronic matrix element with different 
Lorentz structures) whereas at zero recoil one must evaluate all these corrections to a 
single form factor. Furthermore, from an experimental point of view the non-recoil point 
is also more favourable. For example, a possible choice might be the maximum recoil 
point of the J/ip at w — 1.26 (q 2 ~ 0), where several models predict the wave function 
overlap f37|. Nevertheless, notice from figure 3 that, despite smaller vertical (statistical) 
error bars, horizontal (systematic) ones are larger, hence limiting a close access to this 
point and thereby making more advantageous the extrapolation to the w — 1 end-point. 

Lastly we make some final comments: 

• From the experimental viewpoint, observe that if the secondary vertex resolution 
could be improved and/or the lifetime of the B c was closer to 1 ps than to 0.5 ps, 
the experimental accuracy near zero recoil improves considerably. Thus, we show in 
figure 4 possible plots for different angular resolutions. 

• It is important to stress that the quantity to be experimentally measured with a large 
accuracy from B c —>J/i() fi + i> decays can be written as the linear combination: 

r 1/2 x t/ 12 (1) x\V cb \ 



In summary, we suggest to keep an open mind on the possibility of an alternative 
determination of \V c b\ from semileptonic decays of B c mesons at LHC. Even if at the time 
LHC will start to run, B factories would have already provided a (still) more precise value 
of \Vcb\ than at present via the semileptonic B decay, an independent determination of it 
should bring a valuable cross-check. On the other hand, a lot of activity is being devoted 
to the analysis of inclusive production of prompt charmonium resonances at Fermilab. 
Therefore we conclude that if the fragmentation approach to describe the production 
of heavy quarkonia at the large transverse momentum domain is accurately tested (and 
tuned) from experimental data, B c semileptonic decays could be competitive with the B 
ones. 

Alternatively, one can turn the question round and consider \V c b\ as a well-known pa- 
rameter thus verifying QCD calculations, either at the production or at the weak decay 
level. Indeed, fragmentation into heavy quarkonia offers an interesting check of perturba- 
tive QCD and a deep insight into the nonperturbative dynamics in the hadronic formation. 
Besides, a precise knowledge of the B c production rate is a necessary condition for the 
experimental measurement of the absolute branching fraction of any of its decay modes. 
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/V„„/f(w) 




Figure 1: V^b| distribution derived from the (/ 2 -spectrum of the decay B^D*lv measured by ARGUS 

|| and r} A £{w)-\Vcb\ measured by CLEO (T|. 



16 



q 2 (GeV 2 ) 




20 25 30 35 S 23 ( Gey2 ) 



Figure 2: Dalitz plot corresponding to the decay B c ^J/ipfi + u. The leptonic momentum squared is given 
by q 2 = {Pfj, + Pv) 2 = (pb c — Pj/ip) 2 , and £23 = (j>j/$ +Pn) 2 - The horizontal lines correspond to constant 
w = VB c -vj/ip. The non-recoil point of the J/ip in the decaying rest frame has coordinates S23 — rnB c Tnj/^ ~ 
19.5 GeV 2 and q 2 = {piB c —fnj/^) 2 — 10.1 GeV 2 . Notice the linearity between the velocity transfer w and 
q 2 as a kinematic variable characterizing the form factors variation for the hadronic transition. 
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Figure 3: Hypothetical r/i2(w)-\V c b\ distribution of the decay B C -^>J j^^v for a 500 events sample tentatively 
assigned to an integrated luminosity of 20 (i.e. two years of running at low luminosity). The available w 
range extends over the interval (1, 1.26). It was assumed an spatial resolution for flight-path reconstruction 
of a ~ 100 /im and a cut of n = 7 lifetime units (59q ~ 56 mrad). Horizontal bars show systematic errors 
basically limiting the number of experimental points. Vertical bars take into account both the statistical 
fluctuations of the sampled events in the bin (distributed according to the expected differential rate for a 
pseudoscalar to vector semileptonic decay ||]), and the uncertainty coming from contamination of higher 
resonances in the B c decay. Points line up with r]' 12 se t equal to —2. 
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Figure 4: The same as in figure 3 but with different values for 59q corresponding to variations of a between 
50 and 100 fim and r between 0.5 and 1 ps. The improvement on the number of points near zero recoil can 
be considerable. 
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